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A family of non-conjugate chaotic maps generalizing the 
well-known logistic function is defined, and some of its basic 
properties studied. A simple formula for the Lyapunov expo- 
nents of all the maps contained in this family is given based on 
the construction of conjugacies. Moreover, it is shown that, 
despite the dissimilarity of their polynomial expressions, all 
the maps possess the same invariant density. Other algebraic 
properties of the family, which shows some relationship with 
the set of Tschebysheff polynomials, are also investigated. 

PACS numbers: 05.45.-a, 05.45.Ac 



Simple iterated maps, such as the Baker map and 
the logistic map, are the subject of constant fas- 
cination. Part of the interest for these systems is 
linked to the fact that they provide an easy and 
pedagogical way to understand how complex and 
chaotic behavior can arise from simple dynami- 
cal models. Even more remarkable, yet, is the 
fact that studies of low-dimensional maps have 
proven to be fruitful in understanding the ba- 
sic mechanisms responsible for the appearance of 
chaos in a large class of dynamical systems (e.g., 
differential flows, high-dimensional maps). One 
paradigm example of such mechanisms is the so- 
called period-doubling cascades of fixed points, 
encountered qualitatively in many physical sys- 
tems of interest. In this paper, we enlarge the 
set of maps known to be chaotic by presenting a 
generalization of the logistic map. The general- 
ization, more precisely, enables us to construct an 
infinite number of one-dimensional maps which 
are chaotic in the sense that they all have posi- 
tive Lyapunov exponents, and possess at least one 
orbit that is not asymptotically periodic. 



I. INTRODUCTION 

In the theory of nonlinear systems, the logistic map 

Xn+l= fr{Xn)^rXn{l~Xn), (1) 

with < r < 4, e [0, 1], and n = 0, 1, 2, . . ., is well- 
known to provide one of the simplest example of what is 
referred to as a chaotic system. That is, a system which, 
for some range of its parameters, possesses at least one 



bounded orbit {xq,xi, . . .} such that (i) {xq,xi, . . .} is 
not asymptotically periodic and (ii) the Lyapunov ex- 
ponent, defined in the context of Eq.(P by the usual 
relation 



N-l 



\{xo) 



— lim 



(2) 



is greater than zero [y^]. For the logistic map at r — A, 
specifically, X{xo) is positive and equals In 2 for almost all 
xq. This value can be easily calculated by noting that the 
logistic map is conjugate to the tent map or by having 
recourse to the fact that the dynamical system defined 
by Eq. (|l|) is ergodic for this particular value of r . In 
this latter case, one is justified to calculate analytically 
the Lyapunov exponent as an ensemble average 



A= / p{x)lii\dxfr{x)\dx, 
Jo 



the function 
p{x) 



7Ty^x{l — x) 



x e [0, 1] 



(3) 



(4) 



being the invariant density of the map satisfying p{A) = 
p{f~^{A)), where f~^{A) is the preimage of an arbitrary 
subset A of the unit interval |^J|] . 

In this paper, we shall focus on another property of 
the logistic map at r = 4, namely the existence of the 
following closed-form solution 



Xn = Sin2(2"7r0o), 



(5) 



where 9q — ir ^ arcsin(Y^^). The importance of this for- 
mula relies evidently on the fact that be evaluated 
directly, for any initial point xq, without actually com- 
puting the intermediate values xi,X2, ■ ■ ■ ,Xn-i- In the 
remaining of this work, we follow the steps used in the 
derivation of the above formula to construct similar ex- 
pressions involving general terms of the form sit?{N^^6). 
In doing so, we shall see that a whole new family of 
polynomial maps can be defined which, essentially, gen- 
eralize the logistic map by preserving the density Eq. (H) , 
and whose Lyapunov exponents can be calculated easily. 
Other properties of these newly defined dynamical maps 
are also studied. In particular, we shall point out some 
similarities between the maps and the Tschebysheff poly- 
nomials, and prove, finally, that the family of maps as a 
whole is closed under the composition of functions. 



1 



II. DEFINITION OF THE LOGISTIC FAMILY 

The explicit solution of the logistic map at the particu- 
lar value r — A can be derived by substituting in /4(x„) — 
4x„(l — Xn) the change of coordinates Xn = sin^(7r6'„), 
valid for Xn € [0, 1], and by using the trigonometric iden- 
tity 



4sin2 ^(l - sin^ 0) = sm'^{29) 



(6) 



in order to obtain sin^(7r(?„+i) = sin^(27r0„). This is 
equivalent to the map = 20„niodl, which has the 
explicit solution 



9n = 2"6lomodl. 



(7) 



Hence, the complete solution in terms of the coordinate 
X must correspond to Eq.(||). 

Following these steps, one can imagine to construct a 
map having the solution Xn = sin(3"7r0o) by expressing 
sin^(30) in terms of sin^(6') with the identity 



sin^iSe) = 16 sin^ 61-24 sin"* 61 + 9 sin^ 9. 



(8) 



In this case, one gets g{x) = 16x^ — 24a;^ -I- 9a; as a possi- 
ble dynamical map on [0, 1] having the required solution. 
More generally, one can use the recurrence formula 



sin(A^6l) = 2 cos 6* sin[(A^ - 1)^] - sin[(iV - 2)61] (9) 

to define a whole set of maps which express sin^ {N6) in 
terms of sin^ 9. This set, the sine functions set, is defined 
specifically as 

S = {SN{x):N = 1,2,3,...}, (10) 

where SN{x) = [s n {^/x)]'^ , and 

si{x) = X 



S2 



sn 



(x) = 2x\/l - x^ 



(x) = 2\/\ - x^SN-iix) - sn-2{x). 



(11) 



Note that the intermediate functions S]y(x) are expressed 
as such in order to verify Eq. (^) with the variable change 
X = sm9. 



Using these definitions, we calculate the first five func- 
tions of S listed in table 1. Obviously, by construction 
of S, S2{x) is the logistic equation itself with parame- 
ter r = 4. From a more general perspective, it can also 
be seen that the maps SN{x) are degree N polynomi- 
als whose leading coefficient, i.e., the coefficient of the 
highest degree term, is equal to 4^~^ in absolute value. 
These two results, satisfied by any function SN{x), is 
proven more formally in ref. 0|. Note that the latter 
property allows us to extend the similarity with the lo- 
gistic map by parameterizing the functions of iS in the 
following manner 



SNrix) = r 



SNjx) 



(12) 



with < r < 4 . We call the set of functions 
{SNrix) : < r < 4^-i}, the family of SN{x), which 
can be characterized numerically by bifurcation diagrams 
and Lyapunov spectrums such as the ones shown in fig- 
ure 1. 
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FIG. 1. (a)-(b) Bifurcation diagrams of S3{x) and S4:{x) 
(c)-(d) Corresponding Lyapunov spectrums. 



III. PARTICULARITIES OF S 



SN{x) 

Sl{x) = X 

S2{x) = -4x2 + 4a; 

S3{x) = I6x^ - 24x2 + 9x 

S4{x) = -64x* + 128x3 - SOx^ + 16x 

5'5(x) = 256x^ - 640x^ + 560x^ - 200x2 -|- 25x 



Table 1. First five members of S. 



Many of the interesting properties of the logistic map 
at r = 4 can be investigated more intuitively by making 
explicit the fact that Eq.(0) is equivalent to a shift map 
S on the binary expression of ||l|- Indeed, if we express 
9q as a binary number 

^ b 

9o = 0.bohb2... = Y.^' e{0,l}, (13) 

1=0 

then applying Eq. (^ to 6*0 is equivalent to shifting all the 
bits of 9o to the left and dropping the integer part. In 
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other words, 



— O-bnbn+lbn+2 ■ ■ ■ 

= 5"(0.5o5i&2...), 



(14) 



where S*" = 5o S*""^ for n > 1, and 5*^ = 5'. Not surpris- 
ingly, the same is true for the maps SN{x), since Eq.(0) 
was the guidehne in defining the family S. However, in 
the case of SN{x), the shift map to consider takes effect 
on 6*0 written in base N. This follows from the following 
result which generalizes effectively the solution of Eqs.(||) 
and (0). 

Theorem 1. Let {xq, xi,X2, ■ ■ ■} be the orbit of xq un- 
der SN{x). If we write Xn — sin^(7r0„), then we have 



that 



^n+l 



iV"6'o modi. 



(15) 



where, as usual, = tt^^ arcsin(y^). 

We omit the proof of this theorem as it follows directly 
from the next lemma. 

Lemma 1. Consider SN{x) as defined previously. We 
have that 



SN{sva^e) = SIT? [NO). 



(16) 



Proof: The result is obvious for = 1 and N — 2. 
Suppose Eq.(16) true for iV — 1 and N — 2, that is to say 



s^_i (sin e) = S[N -1] (sin^ 9) 
= sm^[{N - 1)6], 

and SN -2{sme) = sm[{N - 2)6]. Then, for 9 G 

SNisiu^ 9) = (sin 61) 



= [2Vl-sin^6lsAr_i(sin6') - SAr_2(sin 6*)]^ 
= [2 cos6lsin((7V - 1)9) - sin((iV - 2)6)f 
= s\T?{N9), 

where we have used the identity (^). □ 
Note that we could have proceeded to a similar general- 
ization of the logistic map using cosine functions instead 
of sine functions, while preserving its shift property. One 
possible way of achieving this is to define the cosine func- 
tions set 



C = {CN{x) : N = 1,2,3, . ..}, 
where CN{x) — [cn{\/x)Y, and 



(17) 



ci(x) = X 
C2{x) = 2x^ 



1 



Cn{x) = 2xCN~l{x) - CN-2{x). 



(18) 



Contrary to the SAr(a;)'s, the functions cn{x) have the 
interesting property that they are polynomials of degree 



N . In fact, the set {cn{x) : N coincides with the set 
of Tschebysheff polynomials on the unit interval 1^ , the 
latter set satisfying the exact same recurrence formula as 
Eq. (|l8|) . We thus have that {cn{x)} must constitute a 
set of orthogonal polynomials, i.e., 



cn{x)cni(x) dx — Sn,n' 



(19) 



for all integers N and N', where Sij is the delta- 
Kronecker function. This fact can be further proved using 
the property Cjv(cos0) — cos{N9), well-known to be sat- 
isfied by the Tschebysheff functions. Note that {sn{x)} 
is also a set of orthogonal functions; its members satisfy 
indeed the relation sn {sin 9) = sin{N9). In the remain- 
ing of this work, we shall restrain our study to the set S, 
since the maps SN{x) are directly related to CN{x) by 
the expression 



SN{x) = 



CN{x), 
1 - CN{x), 



for N odd 
for N even. 



(20) 



Hence, as far as their dynamics are concerned, the func- 
tions SN{x) and CN{x) are totally equivalent. 



IV. CONJUGACIES 

The analysis of the chaoticity properties of a map / is 
greatly simplified by studying conjugate maps of / which 
are obtained by applying a global change of variables. 
Recall that two maps f : I ^ I and g : J J are con- 
jugate if there exists a homeomorphism, i.e., a bijective 
and continuous map H : I J such that H o f = g o H . 
The function H is called a conjugacy. In the context of 
S, a possible conjugate function of SN{x) can be con- 
structed as follows. Let TN{x) be a piecewise linear 
function (a generalized tent map) defined on subinter- 
vals [k/N, {k + 1)/N] of [0, 1] by setting 



TN{x) 



Nx - fc, 
-Nx + k 



1, 



for k even 
for fc odd, 



(21) 



with fc = 0, 1, . . . , AT - 1. 

Theorem 2. SN{x) is conjugate to TN{x) with conju- 
gacy H{x) = sin^(7ra;/2). 

Proof. First note that H{x) is both continuous and 
bijective on the interval [0, 1]. Now, on the first hand we 
have that SN{H{x)) = SN{sm^{TTx/2)) = sin^{NTTx/2). 
On the other hand. 



H{TN{x)) 



s\v?[Nttx/2 — k-K 12), for k even 

sin2(-iV7r2;/2 + (fc-M)7r/2), for fc odd 

(±sin(A^7ra;/2))^ 

sin2(iV7rx/2). (22) 



Thus, we have proved that SN o H ^ H o TN for all 
integers N. □ 
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Figure 2 depicts the graphs of SN{x) and the corre- 
sponding TN{x) for = 3,4. {T2{x) is only the tent 
map since, as we mentioned, 6*2(2;) is the fogistic map.) 
In each case, note that SN{x) is a "smooth" version of 
TN{x). 

Remark. The fmictions CN(x) are pairwise conjugate 
to the functions TN{x) with conjugacy map H{x) — 
cos^(7ra;/2). Evidently, since the conjugacy is an equiva- 
lence relation, this has the consequence that any CN{x) 
is conjugate to SN{x). 
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FIG. 2. Graphs of the generalized tent maps r3(a;) and 
T'iix) which are conjugate to S?>{x) and S'i{x). 



V. ERGODIC PROPERTIES OF S 

In this section, we evaluate the invariant densities of 
the maps SN{x) and their Lyapunov exponents based 
on the conjugates found in the previous section. As one 
would expect from the similarity of the graphs of T'iix) 
and r4(x) in figure 2, the functions TN{x) should exhibit 
a unique and constant invariant density 



Ptn{x) = 1, 



X e [0, 1] 



(23) 



just as it is the case for the tent map. This can be proven 
using Adler and Bowen's results on Markov transforma- 
tions H], (see also Q). At this point, the invariant den- 
sity of the functions SN{x) can be found using the trans- 
formation formula 



PS Nix) = PTNiy) 



dy 
dx 



where y = § arcsin(-yx), so as to obtain 

1 

PS Nix) 



7r-y/x(l — x) 



(24) 



(25) 



The unicity of psNix) is assured by the fact that pTNix) 
is also unique, which means that both functions SNix) 
and TNix) must be ergodic |,|,||. 



Now, to evaluate the Lyapunov exponent of SNix) we 
may use the fact that \dxTNix)\ = N for almost all x G 
[0, 1] to infer that A(a;o) = IniV almost everywhere in the 
case of TNix). Accordingly, since Lyapunov exponents 
are invariant under smooth and differentiable coordinate 
transformations [|l0|, wc have the following theorem. (A 
more extensive proof of this result, which takes care of 
the pathological points where d^TNix) is not defined, is 
contained in ref. 0.) 

Theorem 3. The Lyapunov exponent of SNix) is IniV 
almost everywhere (with respect to the invariant measure 

PSNix)). 

The above theorem shows that the members of S are 
non-conjugate to each other simply because they pos- 
sess different Lyapunov exponents. It also shows that 
iS\{51}, and consequently the set of Tschebysheff poly- 
nomials, are sets of chaotic maps. Indeed, A = IniV > 
for TV > 1, and by using the shift property of SNix) we 
can choose xq = sin^(7r^?o): with 9q irrational, to build an 
orbit that is not asymptotically periodic. Another way 
to convince ourselves that all the polynomials in S have 
chaotic orbits is to use the celebrated result "period-3 im- 
plies chaos" [0 , and find an initial point of period 3 for 
each SNix). For instance, for a iV > 1 let xq = sin^(7r0o) 
where 



1 



0.001001.. 



(in base N). (26) 



" iV3- 1 

Again, using the shift map property, we must have 



Xo 


= sin^(7r 


0.001001 . 


• •) 


Xl 


= sin^(7r 


0.010010. 


■ •) 


X2 


= sm (tt 


0.100100. 


• •) 


X3 


= sm (tt 


0.001001. 


..) ^ xo 



(27) 



We thus extended the chaoticity properties of the logistic 
map to an infinite family of polynomials. 



VI. ALGEBRAIC PROPERTIES 

To complete the study of the properties of S, we now 
deduce that it is an abelian monoid with respect to the 
composition of functions (o). A monoid, precisely, is a 
non-empty set M together with a binary associative op- 
eration, say *, such that x * y £ M for x,y G M. There 
must also be an element e € M, called the identity ele- 
ment, for which x^e = e^x = x for all a; € M . Moreover, 
a monoid is called abelian if the binary operation is com- 
mutative | [Tl[ |. In S, the identity element is 51 (x) = x. 
Also, the composition of function is clearly associative. 
Now, to prove that iS is indeed abelian monoid, we verify 
that it is closed under composition and that this com- 
position is commutative, a condition that is not verified 
in the case of composition of general functions. How- 
ever, before we do so, we present next a new expression 
of SNix) on the unit interval. 
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Lemma 2. For all SN{x) and x e [0, 1], 
SN{x) = sin^ (A^ arcsin v^). 

Proof: Let x G [0, 1]. There exists a. 9 e [0,7r/2] such 
that X — sin^ 0, and thus 9 = arcsin y/x. Now, from 
Lemma 1 we have SN{x) = S'7V(sin^ 6*) = sin^(7Ve') = 
sin^ {N arcsin ^/x) . □ 

Theorem 4- (Monoid property) Let iVi and be any 
positive integers. We have that 



SNi o SN2 = SN2 o SNi = S[NiN2]. 



(28) 



Proof: For A^i and N2 given, consider SNi{x) and 
SN2{x). Then, for any x e [0,1], we obtain from 
Lemma 2 



SNi{SN2{x)) = sin^lNiaTCsm{y/SN2{x))] 

= sin^[A^i arcsin(sin(A''2 arcsin -^/a;))] 



■(iViiV2 arcsin ^/x) 



= S[NiN2]ix) 



(29) 



Obviously, S[NiN2]ix) = S[N2Ni]{x), so the composi- 
tion is commutative. □ 
As a direct consequence of the monoid property, k- 
periodic points of a certain polynomial SN{x) can be 
looked at as fixed points of the function SM{x) where 
M = . Furthermore, a polynomial SN{x) of very 
high degree can be computed easily by decomposing its 
expression using lower degree polynomial of the family 
S. Exphcitly, consider SN{x) £ S. We say that SN{x) 
is a prime element of 5 if A'^ is a prime number. Us- 
ing this definition, we have as a result of Theorem 4 and 
the Fundamental Theorem of Arithmetic that any poly- 
nomial SN{x) must possess a unique decomposition in 
prime elements of S. 



VII. FINAL REMARKS 

To conclude, note that our study of the sine func- 
tions, written in the form SN(x) = sin^ (A^ arcsin y^) , 
have been restricted to positive integers N. In a simi- 
lar manner, it could be interesting to investigate func- 
tions of the type Sa{x) = sin^ (a arcsin -yx) with a 
real. One observation about this extra generalization 
is that, as for Sl{x), the function Sa(x) does not ex- 
hibit chaotic properties for < a < 1. The function 
S^{x) = sin^ ( i arcsin Y^) , for example, is conjugate to 
g{x) — x/2, and has all of its orbits attracted to a: = 0. 
Yet, this is not surprising since the Lyapunov exponent 
of this map must be ln(l/2) < 0. This brings us to con- 
jecture that Sa{x) must admit chaotic behavior if and 
only if |a| > 1, considering that the Lyapunov exponents 
of Sa{x) should be In a. A complete proof of this result, 
however, cannot be given here using the same symbolic 
dynamic approach used for SN{x), for the simple reason 
that the expression of a point in "base TV" makes sense 
only if N is an integer greater than 1. 
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